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ABSTRACT 

Aims. We examine the effects of the horizontal turbulence in differentially rotating stars on the GSF instability and 
apply our results to pre-supernova models. 

Methods. We derive the expression for the GSF instability with account of the thermal transport and smoothing of the 
/^-gradient by the horizontal turbulence. We apply the new expressions in numerical models of a 20 Mq star. 
Results. We show that if Nq < the Rayleigh-Taylor instability cannot be killed by the stabilizing thermal and fi- 
gradients, so that the GSF instability is always there and we derive the corresponding diffusion coefficient. The GSF 
instability grows towards the very latest stages of stellar evolution. Close to the deep convective zones in pre-supernova 
stages, the transport coefficient of elements and angular momentum by the GSF instability can very locally be larger 
than the shear instability and even as large as the thermal diffusivity. However the zones over which the GSF instability 
is acting are extremely narrow and there is not enough time left before the supernova explosion for a significant mixing 
to occur. Thus, even when the inhibiting effects of the /^-gradient are reduced by the horizontal turbulence, the GSF 
instability remains insignificant for the evolution. 

Conclusions. We conclude that the GSF instability in pre-supernova stages cannot be held responsible for the relatively 
low rotation rate of pulsars compared to the predictions of rotating star models. 
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1. Introduction 

The comparison of the observed rotation rate of pulsars 
and stellar models in the pre-supernova stages indicate 
that most stars are losing more angular momentum than 
currently predicted QHeger et al. 2000} IHirschi et al. 2 004) . 
Normally, the conservation of the central angular momen- 
tum of a presupernova model would lead to a neutron 
star spinning with a period of 0.1 ms, which is about 
two orders of magnitude faster than the estimate for 
the most rapid pulsars at birth. The question has arisen 
whether some rotational instabilities may play a role in 
dissipating the angular momentum. We can think in par- 
ticular of the Golreich-Schubert-Fricke (GSF) instability 
(|Goldreich fe Schubert 19671 IFricke 1968]) . which has a neg- 
ligible effect in the Main-Sequence phase and which may 
play some role in the He-burning and more advanced 
phases (Heger et al. 2000), in particular when there is a 
very steep fi-gradient at the edge of the central dense core. 
This instability is generally not accounted for in stellar 
modeling. The aim of this article is to examine whether the 
GSF instability is important in the pre-supernova stages, 
when account is given to the effect of the horizontal turbu- 
lence in rotating stars which reduces the stabilizing effects 
of the /^-gradient. 

Sect. 2 recalls the basic properties of the GSF instability, 
Sect. 3 those of the horizontal turbulence. The effects of 



turbulence on the GSF instability are examined in Sect. 4. 
Sect. 5 show the results of the numerical models. Sect. 6 
gives the conclusion. 



2. The GSF Instability and Solberg-Hoiland 
Criterion 

2.1. Recall of basics 

A rotating star with a distribution of the specific angu- 
lar momentum j decreasing outwards is subject to the 
Rayleigh-Taylor instability: an upward displaced fluid el- 
ement will have a higher j than the ambient medium and 
thus it will continue to move outwards. In radiative sta- 
ble media, the density stratification has a stabilizing effect, 
which may counterbalance the instability resulting from 
the outwards decrease of j. In this respect, the /^-gradient 
resulting from nuclear evolution has a strong stabilizing 
effect. The stability condition is usually expressed by the 
Solberg-Hoiland criterion, given in the first part of Eq. ([!]) . 

The GSF instability occurs when the heat diffusion by 
the fluid elements reduces the stabilizing effect of the en- 
tropy stratification in the radiative layers. The account of a 
finite viscosity v together with thermal diffusivity K influ- 
ences the instability criteria (iFricke 19681 [Xcheson 1978). 
These authors found instability for each of the two condi- 
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where Nj, ad is the adiabatic thermal term of the Brunt- 
Vaisala, (BV) frequency and the rotational contribution 
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The viscosity v = (l/3)v£ represents any source of viscosity, 
including turbulence, w is the distance to the rotation axis 
and z the vertical coordinate parallel to the rotation axis. 
The thermal diffusivity K is 
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where the various quantities have their usual meaning. 



— The first inequality in Eq. ([!]) corresponds to the convec- 
tive instability predicted by the Solberg-Hoiland crite- 
rion with account for the efficiency factor T — v£/(6K) 
which takes into account the radiative losses. For Nq < 
0, a displaced fluid element experiences a centrifugal 
force larger than in the surrounding and further moves 
away. The first criterion in Eq. (fTl) expresses that insta- 
bility arises if the T gradient, with account for thermal 
and viscous diffusivities, is insufficient to compensate for 
the growth of the centrifugal force during an arbitrary 
small displacement. 

— The second inequality in Eq. ([I]) expresses a baroclinic 
instability related to the differential rotation in the di- 
rection z. If a fluid element is displaced over a length Sz 
in the z direction, so that dil/dz-Sz > 0, the angular ve- 
locity of the fluid element is larger than the local angular 
velocity. The excess of centrifugal force on this element 
leads to a further displacement and thus to instability. 
It has often been concluded from this second criterion 
that only cylindrical rotation laws are stable (solid body 
rotation being a peculiar case). This is not correct, since 
viscosity is never zero. In particular the horizontal tur- 
bulence produces a strong horizontal viscous coupling, 
with a large ratio vjK, which does not favor the insta- 
bility due to the second condition in Eq. (JfJ. 

Numerical simulations of the GSF instability 
( jKorycansky 1991] ) show that the GSF instability de- 
velops in the form of a finger-like vortex in the radial 
direction, with a growth rate comparable to that of the 
linear theory. 



2.2. The fi gradient and the GSF Instability 

In the course of evolution, a fi gradient develops around 
the convective core (there the gradients are also large). 
The \i gradient produces stabilizing effects. Endal and Sofia 
(|1978|) in their developments surprisingly use the same de- 
pendence on the /i-gradient as for the meridional circula- 
tion (see also |Heger et al. 2 000 ) . They apply a velocity of 
the GSF instability in the equatorial plane given by 



^gsf = 2 — — — U 2 {r) 
Hi dmr 



(4) 



where U 2 (r) is the radial component of the velocity of 
meridional circulation and Ht and Hj are respectively the 
scale heights of the distributions of T and specific angular 
momentum j. 

Let us focus on the first criterion i n Eq. (jlj), it becomes 
in this case ( |Knobloch fc Spruit 19831 ITalon 1997|) 
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(2) is the particle diffusivity, either molecular or radiative. 

It is generally of the same order as the viscosity v, thus the 
stabilizing effect of the \i gradient is not much reduced by 
the diffusion of particles. Thus, when there is a significant 
fx gradient, it generally dominates and tend to stabilize the 
medium. This is why the GSF instability is generally of only 
limited importance in regions with Nq < surrounding the 
stellar cores in advanced phases. The occurrence of horizon- 
tal turbulence, however, greatly changes the above picture, 
because it is anisotropic and produces a very large particle 
diffusivity, thus reducing the effect of the jji gradient. 



3. The Coefficient of Horizontal Turbulence in 
Differentially Rotating Stars 

The importance of the horizontal turbulence in differen- 
tially rotating stars was emphasized by Zahn (1992). There 
are a number of observational effects supporting its ex- 
istence, in particular the thinness of the solar tachoclinc 
( |Spiegel fc Zahn T9 92 ) , the different efficiencies of the 
transport of chemical elements and of angular momen- 
tum as well the observations of the Li abundances in solar 
type stars ( |Chaboyer et al. 1995a[ [Chaboyer et al. 1995b[ ). 
In massive stars, the horizontal turbulence increases the 
mixing of CNO elements in a favorable way with respect to 
observations (IMaeder 2003|) . 

A first estimate of the coefficient of horizontal tur- 
bulence was proposed by Zahn (1992). A second better esti- 
mate was based on laboratory experiments with a Couette- 
Taylor cylinder. It gives in a differentially rotating medium 
(jRichard fc Zahn 1999|IMathis et al. 2004p . 
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The latitudinal variations of the angular velocity are 
of the form fi(r,i?) = H(r) + fi(r,tf) = O(r) + 
na(r) (P 2 (??) + 5).nisthe average on an isobar, while O2 
expresses the horizontal differential rotation (Zahn 1 1992 . 
IMathis fc Zahn 2004p . 
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with a = h d M r n ) . In a star with shellular rotation, one 

_2_ a In r ' 

has f^2 <C 0,(r). The diffusion coefficient of horizontal tur- 
bulence (which is also the viscosity coefficient) becomes 
([Mathis et al. 2004|) . 



D h = u h = -pr 2 \n 2 \ 



1/2 



1/2 [r\2V 2 -aU 2 \} 1/2 



(8) 



Raphael Hirschi 1 ' 2 , Andre Maeder 3 : The GSF instability and turbulence 



3 



where U 2 and V% are the vertical and horizontal components 
of the velocity of meridional circulation and a is the same 
numerical factor as in Eq. (Fn). 

The above diffusion coefficient (Eq. [8]) derived from lab- 
oratory experiments is essentially the definition of the vis- 
cosity or diffusion coefficient, if the characteristic timescale 
of the process is equal to l/(/3 2 ), i.e. 
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with ^(jiff 
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with £ ~ r. This relation implies that only the degree of 
the differential rotation in 1? determines the importance of 
horizontal turbulence. However, the motions on an isobar 
in spherical geometry are not necessarily the same as in the 
Couette-Taylor experiment of rotating cylinders, which is 
only a local approximation of the horizontal shear on a tan- 
gent plane. If the horizontal turbulence is rather related to 
the differential effects of the Coriolis force (Maeder 2003) 
which acts horizontally, i.e. idifi ~ r/(Q 2 V2) 1 / 2 , one ob- 
tains the following coefficient 

v h = Ar\rU{r) V 2 [2V 2 ~ aU 2 ] J with A < 0.1 . (10) 



This expression, despite its difference with respect to Eq. 
{[sj) , leads to similar numerical values for the horizontal tur- 
bulence in stellar models (jMathis et al. 2004]) . while the 
original estimate (Zahn 1992) leads to a coefficient _D n 
smaller by four orders of a magnitude. 

The expression of requires that we know the verti- 
cal and horizontal components U 2 and V 2 of the velocity 
of meridional circulation. If not, some approximations are 
given in the Appendix. 

4. The Horizontal Turbulence and the GSF 
Instability 

We examine what happens to the condition ^ or Solberg- 
Hoiland criterion in case of thermal diffusivity and horizon- 
tal turbulence. For that let us start from the Brunt -Vaisala 
frequency in a rotating star at colatitude d 



The ratio T/(r + 1) in Eq. ( 12 ) is the fraction of the energy 
transported. 

The GSF instability problem is 2D with two different 
coupled geometries: the cylindrical one associated to the 
rotation with the restoring force being along e s and the 
spherical one where the entropy and chemical stratication 
restoring force is along e r that explains the sin d in Eq. (11), 



which gives the radial component of the total restoring 
force. The following formula for Nq in spherical geometry 
in the case of a shellular rotation ft(r) can be obtained: 
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then introducing sp : 
der to simplify the problem, we will focus on the equatorial 
plane (i9 = 7r/2), in which case we simply have: 



starting with Eq. (|2j: - jg ds 
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2ft 



din ft 
d In r 



(15) 



The horizontal turbulence also makes some exchanges 
between a moving fluid clement with composition given by 
/lint and its surroundings with mean molecular weight ji. If 
fn is the amount of fi transported expressed in fraction of 
the external gradient, one has 



p /i,int 



(16) 



One can also write / M = r^/(r^ + 1), where Y ^ is the ratio 
of amount of /x transported to that lost by the fluid element 
on its way. Thus, one has 



/i,int 
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v£ 
6A^ 



(17) 



to be compared to the first part of Eq. 12 If N 2 < 0, the 
medium is unstable, thus the instability condition at the 
equator becomes 



N 2 = N^ + N 2 + Nq sintf = 
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If it is negative, the medium is unstable. V; n t is the internal 
gradient in a displaced fluid element, while V is the gradi- 
ent in the ambient medium. These gradients obey to the 
relations (Mae der 1995|) 
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For a fluid element moving at velocity v over a distance £, 
T = Pe/6 = v£/(6K), where Pe is the Peclct number, i.e. 
the ratio of the thermal to the dynamical timescale. T is the 
ratio of the energy transported to the energy lost on the way 
£. The horizontal turbulence adds its contribution to the ra- 
diative heat transport and T becomes (jTalon fc Zahn 1997ft 
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The situation is similar to the effect of horizontal turbulence 
in the case of the shear instability (Talon & Za hn 1997)) . 

The turbulent eddies with the largest sizes x — v£/6 
are those which give the largest contributio n to the vertical 
transport. For these eddies, the equality in ( 18 ) is satisfied, 
which gives 



K + D h 



^T,ad 



— N 2 







(19) 



The diffusion coefficient by the GSF instability is -Dgsf = 
(l/3)v£ = 2x, obtained from the solution of this second 
order equation, which may also be written, 



K 2 d 



[N 2 d D h + N 2 (K + D h ) + N 2 (K + 2D h )) x + 
N 2 (D h K + D 2 ) = 0. 
We notice several interesting properties. 



(20) 
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$ = 150 kms- 1 , Z M = 0.002 
He-burning 




M = 20 AI V , Vitf = 150 kms" 1 , Z m = 0.002 
Ne-burning 




Fig. 1. Properties of a 20 M model with Z = 0.002 and 
an initial rotation velocity of 150 km s _1 during the He- 
burning phase, when the central He content is Y c — 0.543 
and the actual mass 19.795 M . a) The top left panel il- 
lustrates the various diffusion coefficients as functions of 
the internal mass. The grey areas correspond to convective 
zones, b) The top right panel shows the profile of the angu- 
lar velocity il-gradient (dlnO/dlnr) + 2. A negative value 
of this term means instability, c) The left bottom panel 
shows the various iV 2 . d) The right bottom panel shows 
the profile of fl and its ratio to the local critical angular 
velocity fJ cri t- 



1. If < 0, from Eq. (|19|) we see that the GSF instability 



is present in a radiative medium whatever the \x- and 
T-gradients are. Thus, these gradients cannot kill the 
turbulent transport by the GSF instability. However, 
the size of the effects has to be determined for any given 
conditions. 

2. If the diffusion coefficient -Dqsf by the GSF instability 
is small with respect to K and Z?h, we have 



D, 



GSF 
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( 



\(K+D h ) + AT 



(21) 



The assumptions D GSF <C K and -Dqsf "C Ai ar e 
likely, at least at the beginning of the GSF instability 
when Nq starts becoming negative. Nevertheless, these 
assumptions need to be verified for the cases of interest 
in the advanced stages. 
3. If N 2 ^> Njt ad , as is the case in regions surrounding 



stellar cores, we get from Eq. ( 19 ) 
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Fig. 2. Same as for Fig. 1 during the phase of central 
Ne-burning. The actual mass is 19.412 Mq, the central 
Ne content is X( 20 Ne) =0.261, the central O-content is 
X( 16 O)=0.726. 



No assumption on the size of Dgsf is made here. Due 
to the fast central rotation, and the f2-gradient in 
regions close to the central core may be large, thus pos- 
sibly favoring a significant Dqsf- 

For more general cases, the simple solution of the second 
order equation (f20j) has to be used. Most critical of course 
are the values oflV^ and N 2 , which take large values in a 
narrow region surrounding the central core in the helium 
and more advanced evolutionary stages. 



5. Rotating Stellar Models in the Pre-Supernova 
Stages 

In order to examine quantitatively the importance of the 
GSF instability, we calculate the evolution all the way from 
the Main Sequence to the Si burning stage of a 20 M star 
with an initial rotation velocity of 150 km s _1 with a metal- 
licity Z=0.002 typical of the SMC composition. We make 
the choice of this composition, because the internal f2- 
gradients are steeper at lower Z ( |Maeder k, Meynet 2001] ), 
which would favor the GSF instability. Some data for an- 
other 20 M model with an initial rotation of 300 km s^ 1 
are also given. Equation 
curence of the GSF insta 



20 1 was used to determine the oc- 
Dihty and the value of Z?gsf- The 
above expression ( 10 ) for is used. The nuclear network 
in the advanced phases is the same as in previous models 
(THirschi et al. 2004p . 

Figure [T] shows in four panels the main parameters dur- 
ing the first part of the phase of central He-burning. We 
first notice in panel d) the building of a f2-gradicnt at the 
edge of the convective core with a difference of fl by about 
a factor of 20. This makes dmfi/dlnr + 2 < in most of 
the region between the edge of the convective core at 2.9 
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i = 150 kms- 1 , Z M = 0.002 
O-burning 



= 20 M B , »,„ = 300 kms- 1 , Z m = 0.002 
He-burning 





Fig. 3. Same as for Fig. 1 during the phase of central O- 
burning. The actual mass is no longer changing, the cen- 
tral O-content is AT( 16 O)=0.692, the central Si-content is 
X( 28 Si)=0.153. 



M Q and the convective H-burning shell at 5.3 M Q as shown 
in panel b) . However Nq remains negligible with respect to 
Nj- and N?. In order to understand why, we need to look 



Fig. 4. Same as for Fig. 1 for an initial velocity of 300 km 
s _1 . The star is in the stage of central He-burning with 



back at Eq. (151: Nq = 2fl (2 + f^). The value of n 2 
in the star is too small to allow a significant value of N^. 
This means in fact that the centrifugal force in the deep 
interior is not strong enough to overcome the stabilizing ef- 
fects of iVy, and iV 2 as shown in panel c). The consequence 
as illustrated in panel a) is that £>gsf remains everywhere 
smaller than -D s hear and is thus not significant. We also no- 
tice that Dgsf is always much smalle r th an and K, 
which permits here the approximation (21 1 made above. 



Figure [2] shows the same plots during the stage of cen- 
tral neon burning. We notice an impressive increase of the 
central angular velocity and a very small Q in the envelope, 
with a difference by a factor of 10 s between the two, justi- 
fying the examination of the GSF instability. The are two 
"fi-walls" , the big one at 7.2 M corresponds to the basis of 
the H-rich envelope, the other one at 4.8 M lies at the ba- 
sis of the He-burning shell. The values of (d In il/d In r) + 2 
become much more negative, however over areas of very 
limited extensions. Again, the value of are negligible, 
in particular compared to the big peak of iV 2 at 4.8 M Q . 
The result is that -Dgsf is always smaller than -D s h car , even 
if very locally it can reach about the same value. -Dgsf is 
always at least two or three orders of a magnitude sma ller 
than Dh and K, permitting here the simplification (21). 



Figure [3] shows the situation in the central O-burning 
stage slightly less than a year before the central core col- 
lapse. Two other small steps in Q have appeared near 
the center, due to the successive "onion skins" of the 
pre-supernova model. We notice some new facts. In line 
with what was already seen for neon burning, the term 



Y" r =0.247. The actual mass is 19.681 M, 



0- 



(d In Q/d In r) + 2 becomes negative only in extremely nar- 
row regions where the GSF instability is acting with a dif- 
fusion coefficient Dgsf larger than in the previous evolu- 
tionary stages. Very locally at the upper and/or lower edges 
of intermediate convective zones, Dgsf may even become 
larger than and K reaching values above 10 8 cm 2 s _1 
(there approximation (21) is not valid!). With less than a 
year left before explosion, the distance over which a signif- 
icant spread may occur is about 10 -3 R©. This is not en- 
tirely negligible in the dense central regions, however this 
remains of limited importance, as shown by panels b) and 
d) where we notice that the fi-walls remain unmodified 
despite the locally large Z? G sf- 

We may wonder whether higher initial rotation veloc- 
ities lead to different results. Figure [2] shows the various 
panels for a similar star in the He-burning phase with an 
initial rotation velocity of 300 km s -1 . We see that the cen- 
tral rotation velocity is about the same as for the previous 
case of lower rotation and, in this stage which determines 
the further evolution, there is no significant difference in 
the various properties. 



6. Conclusions 

We have examined the effects of the horizontal turbulence 
on the GSF instability. This instability is present as soon 
as Nq is smaller than zero, whatever the effects of the sta- 
bilizing /i-gradicnts. 

On the whole, the numerical models of rotating stars 
show that the diffusion coefficient by the GSF instability 
grows towards the very latest stages of stellar evolution, 
however the zones over which it is acting arc extremely 
narrow and there is not enough time left before the su- 
pernova explosion for a significant mixing to occur. Thus, 
even when the inhibiting effect of the /i-gradient is reduced 
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by horizontal turbulence, the GSF instability is unable to 
smooth the steep fi-gradients and to significantly transport 
matter. 

We conclude that the amplitude and spatial ex- 
tension of the GSF instability makes it unable to 
reduce the angular momentum of the stellar cores 
in the pre-supernova stages by two orders of mag- 
nitude. Therefore, other mechanisms such as mag- 
netic fields (jSpruit 2002], |Maeder fe Meynet 2004 



IMathis fe Zahn 2005] IZahn et al. 2007P and gravity 
waves (|Talon fc Charbonnel 20051 IMathis et al. 2008|) 
must be further investigated. 

Appendix: some approximations for meridional 
circulation 



Maeder, A 
Maeder, A 
Mathis, S., 
Mathis, S., 
Mathis, S., 
Mathis, S., 

251, 
Richard 



Meynet, G. (2001) A&A 373, 575 
Meynet, G. (2004) A&A 422, 225 
i., Palacios, A., Zahn, J.-P. (2004) A&A 425, 243 
3., Zahn, J.-P. (2004) A&A 425, 229 
3., Zahn, J.-P. (2005) A&A 440, 653 

i., Talon, S., Pantillon, F.-P., Zahn, J.-P. (2008) Solar Physics 
101 

D., Zahn, J.-P. (1992) A&A 347, 734 
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Spiegel, E., Zahn, J. P. (1992) A&A 265, 
Spruit, H. C. (2002) A&A 381, 923 
Talon, S. (1997) Thesis, Univ. Paris VII, 187 p. 
Talon, S., Zahn, J.-P. (1997) A&A 317, 749 
Talon, S., Charbonnel, C. (2005) A&A 440, 981 
Zahn, J.P. (1992) A&A 265, 115 
Zahn, J.P., Brun, A. S., Mathis, S. (2007) A&A 474, 
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The coefficient -Dqsf requires, because of the horizon- 
tal turbulence, the knowledge of the components U 2 and 
V 2 of the meridional circulation. If the solutions of the 4 th 
order system of equations governing meridional circulation 
are not available, some approximations may be considered. 
We note that the same problem would occur for Eq. ( 4]) by 
Endal and Sofia (1978). As shown by stellar models, the or- 
ders of magnitude of U 2 and V 2 are the same. The numerical 
models give in general V 2 ~ U 2 /3 and \2V 2 — ctU 2 \ ~ V 2 . 
Using these orders of magnitude in Eq. (8j), we get 
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1/2 



rU 2 
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(24) 



For U 2l various expressions can be used taking into account 
the amount of differential rotation ( Maeder 2009 ) . We can 
also get an order of magnitude using the approximation for 
a mixture of perfect gas and radiation with a local angular 
velocity fi(r), ignoring the effects of differential rotation on 
the circulation velocity and the Gratton-Opik term which 
is large only in the outer layers, 
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where the various quantities have their usual meaning. 
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